Spin-triplet f-wave pairing due to three-site cyclic-exchange ferromagnetic interactions 
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Ferromagnetiam and superconductivity in a two-dimensional triangular-lattice Hubbard model 
are studied using the density-matrix renormalization group method. We propose a mechanism of 
the /-wave spin-triplet pairing derived from the three-site cyclic-exchange ferromagnetic interactions. 
We point out that a triangular network of hopping integrals, which is required for the three-site 
cyclic hopping processes, is contained in the (possibly) spin-triplet superconducting systems, such 
as Bechgaard salts (TMTSF)2X, cobalt oxide Nao. 35 C0O2 -1.31120, and layered perovskite Sr2Ru04. 

PACS numbers: 71.10.Fd, 74.20.-z, 74.20.Rp, 75.50.Cc 



Symmetry of Cooper pairs give an invaluable infor- 
mation of our elucidating a superconducting mechanism. 
Both spin and angular momentum are good quantum 
numbers for describing the Cooper pairs. In the spin- 
triplet channel, the total spin is one and the pairing sym- 
metry is of p- or /-wave. The long-standing quest for a 
new example of spin-triplet superconductivity has been 
pursued since the spin-triplet superfluid ^He was discov- 
ered [l[. Recently, three experimental candidates have 
been hand-running recognized to be likely spin-triplet 
superconductivity (see below). They provide an excel- 
lent chance for us to have productive insights into the 
nature of spin-triplet pairing, which is less understood 
than that of the spin-singlet pairing. Spin-triplet super- 
conductivity is thus one of the hottest topics in the field 
of strongly-correlated electron systems. 

From the theoretical aspect, spin-triplet Cooper pairs 
are considered to be formed simply by a ferromagnetic in- 
teraction . The realization of such pairs in a system of 
electrons interacting via spin-independent Coulomb in- 
teraction, namely the Hubbard model, is still a signif- 
icant challenge. The Nagaoka-Thouless 0, and flat- 
band mechanisms are well known as the origin of fer- 
romagnctism. However, both of them are unsuitable for 
explaining the spin-triplet superconductivity since they 
lead not to Cooper pairs but to saturated magnetiza- 
tion. Then one may turn his eyes to the three-site cyclic- 
exchange ferromagnetic interaction which is only the re- 
maining mechanism in the simple Hubbard model @ . It 
brings a couple of electrons to be locally ferromagnetic, 
so that applying the mechanism to the spin-triplet su- 
perconductivity seems to be a natural extension. In this 
letter, we propose a mechanism of the spin-triplet super- 
conductivity derived from the three-site cyclic-exchange 
ferromagnetic interaction. 

We are also highly motivated by the fact that all of 
the three experimental candidates contain a triangular 
network of hopping integrals which is required for the 
three-site cyclic-exchange mechanism. One of the can- 



didates is Bechgaard sahs (TMTSF)2X 0. The crys- 
tal structure consists of well-separated sheets containing 
one-dimensional (ID) TMTSF stacks along the a-axis. 
The sheets are in the a6-plane and the hopping integrals 
along the 6-axis arc about 10 — 20% of those along the 
a-axis Q . The unique structure of the hopping integrals 
can be regarded as an anisotropic triangular lattice. In 
experiment [i| , it has been widely concluded that the su- 
perconducting state is of the spin-triplet pairing: the spin 
susceptibility, from the NMR Knight shift measurement, 
remains unchanged through the superconducting transi- 
tion (Tc) and the upper critical magnetic field exceeds 
substantially the Pauli limit in both a and h directions. 
Also, the enhancement of the spin-lattice relaxation ra- 
tio just below Tc, indicating the presence of line nodes in 
the superconducting order parameter, was found [lo| . In 
theory llj , we seem to have reached a consensus on the 
spin-triplet pairing; nevertheless, the pairing symmetry 
is still open issue [l3, 0, 0, 0, [ll, [13, [3 ■ 



Another experimental candidate is cobalt oxide 
Nao.35CoO2-l.3H2O, which shows a superconducting 
transition at Tc « 5 K [l^ . The conductive ah planes 
consist of edge-sharing CoOe octahedra and each plane 
is strongly separated by Na"*" ions and H2O molecules 
along c axis. The Co ions form a triangular lattice, and 
the system may be regarded as a two-dimensional (2D) 
triangular lattice doped with 35% electrons. Some ex- 
periments have suggested the possibility of an uncon- 
ventional superconductivity with line nodes in the gap 
function [20l. [2ll. [23 . [23j . In particular, spin-triplet p- 
or /-wave pairing state has been argued based on an in- 
variant bahavior of the a6-plane Knight shift across Tc 
and a T^-dependence of 1/Ti below Tc [13]. On the 
other hand, a suppression of Knight shift across Tc for 
both magnetic field parallel and perpendicular to the 
afe-plane suggesting a spin-singlet pairing state has also 
been reported 24|. It is supported by an angle-resolved 
photoemission study, which have observed only hole-like 
Fermi surface aroimd F point [2^. As described above. 
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experimental decision of the pairing symmetry is still 
controversial. Similarly, theoretical studies are not in 
agreement: the d+id- (or p-)wave spin-singlet pairing 
with a triangular-lattice t—J model [2^, 27, 2^, 2^ and 
the /-wave spin-triplet pairing with extended Hubbard 



model [30|, |3l|, |33l have been proposed. Also, the p-wave 
spin-triplet has been suggested phenomenologically [s^l- 
The third candidate is Sr2Ru04, which is a tetragonal, 
layered perovskite system of stacking Ru02-planes. No 
sooner was discovered the superconducting transition at 
Tc ~ 1.5 K [34, the possibility of spin-triplet pairing was 
pointed out [35[. After that, the spin-triplet pairing state 
with d-vector perpendicular to the conducting plane has 
been confirmed by a NMR measurement [s^ . Like Cu in 
the high-Tc superconductors, Ru forms a square lattice. 
If a single-band description of Ru02-plane for so-called 
7 band could be adequate, the system is described as 
a 2D Hubbard model with next-nearest-neighbor hop- 
ping. The next-nearest-neighbor hopping has been esti- 
mated to be 0.3 — 0.4 in units of the nearest-neighbor hop- 
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so that Ru indeed forms a triangular net- 
Based on the quantum oscillation measurement, 
the 7 Fermi surface sheet is a large electron-like cylinder 
with the electron filling n ~ 2/3. At present, there exist 
several theoretical studies where a competition between 
d-wave spin- sing let and /-wave spin-triplet pairing states 
is discussed H SO, El 



In our previous works 43, 44|, we studied the ground- 
state properties of two-chain Hubbard model coupled 
with zigzag bonds. This system is equivalent to a ID 
triangular-lattice Hubbard model. We argued there that 
a cyclic hopping motion of two electrons in each triangle 
yields a ferromagnetic correlation in the strong-coupling 
regime; as a consequence, a spin-triplet superconducting 
state becomes dominant if the following conditions are 
satisfied: (i) The product of three hopping integrals in 
each triangle is positive, (ii) The zigzag bond is rather 
weaker than the intra-chain bond, (iii) The filling is sub- 
stantially away from zero and half fillings. The three-site 
cyclic-hopping mechanism should work even in 2D sys- 
tem, and it is particularly worth noting that the above- 
mentioned three materials completely meet the condi- 
tions. Wc therefore consider a 2D anisotropic triangular- 
lattice Hubbard model. We employ the density-matrix 
renormalization group (DMRG) method [i^ to calculate 
the spin-spin correlation and spin-triplet pair-correlation 
functions. It will thereby be confirmed that the spin- 
triplet pairing occurs predominantly in the /-wave chan- 
nel for a simple 2D Hubbard model. 

The Hamiltonian of the anisotropic triangular-lattice 
Hubbard model is given by 

H=Y1 + H.C.) + uY,n^^rl,l (1) 

where cj^ (cia) creates (annihilates) an electron with spin 
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FIG. 1: (Color online) Schematic representation of the 
anisotropic triangular-lattice Hubbard model. Lattice drawn 
with bold lines is extracted for our calculation of the spin- 
triplet pair-correlation functions. The f-wave pair function is 
also shown. 



a at site i, rii 



is a number operator, and ty is 
hopping integral between sites i and /. The sum {ij) 
runs over nearest-neighbor pairs. We here include two 
kinds of nearest-neighbor hopping integrals t in the x- 
direction and t' otherwise, as shown in Fig. [1] Onsite 
Coulomb repulsion U is set to be \Qt as a typical value in 
the strongly-correlated electrons system. We take i = 1 
as the energy unit hereafter. 

Let us first demonstrate that the cyclic-hopping mech- 
anism indeed leads to ferromagnetic correlation in the 
2D system ([T|). We thus study a finite-size cluster of 
Ix = 6 and ly ~ 8, where the periodic (open) bound- 
ary conditions arc applied for the x- [y-) direction. This 
choice of the boundary conditions enables us to carry 
out sufficiently accurate calculations. We keep up to 
TO = 3400 density-matrix eigenstates in the DMRG pro- 
cedure. In this way, the maximum truncation error, i.e., 
the discarded weight, is less than 1 x 10~^ and the max- 
imum error in the ground-state energy is estimated to 




FIG. 2: (Color online) Average spin-spin correlations <^Sf 
between sites bonded with i', as a function of filling n. 
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We next turn to the long-range behavior and the spa- 
tial symmetry of spin-triplet pair-correlation functions 
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FIG. 3: (Color online) (a) Definition of the rotation angle 6 
and the distance I. For a case that the pair is removed from a 
position i with dashed ellipse, the position of the pair created 
at j -I- Z is denoted by solid ellipse for each 9. (b) log \D{6, l)\ 
with I &t t' = 0.5 for all possible 6. Filled (empty) symbol 
denotes a positive (negative) value of D{9,1). 



be Eo/{lxly) ~ 10 ^. Figure [2] shows the average spin- 
spin correlations l^SfS^) between the neighboring sites 
bonded with t', where (• ■ •) denotes the ground-state ex- 
pectation value. We can see the enhancement of ferro- 
magnetic correlation in a wide range of filling for t' = 0.25 
and 0.5. The suppression of the correlation around 
n — 0.7 implies that the origin of ferromagnetism for 
n > 0.7 is different from that for n ^ 0.7: it comes from 
the cyclic-hopping mechanism around n = 0.5 and from 
the Nagaoka-Thouless mechanism near n = 1. Those 
behaviors are qualitatively consistent with the results in 
the ID triangular-lattice Hubbard model [Si. On the 
other hand, antiferromagnetic correlation is mostly dom- 
inant for t' = 1 since the direct exchange interaction will 
exceed the ferromagnetic ones. We note that the correla- 
tion seems to be slightly ferromagnetic for rt < 0.2, where 
the nearly flat-band system may be realized. So we con- 
firm that in the 2D triangular-lattice Hubbard model the 
ferromagnetic correlation induced by the cyclic-hopping 
mechanism can occur at n 0.5 if the anisotropy is 
rather strong. 
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with a spin-triplet operator — Ci^Ci-^-ri + Ci|Ci_)-rt 
where i + r denotes the neighboring site of i. As in Fig.[T] 
we extract a three-leg ladder from 2D lattice, which is a 
minimal model to generate all possible pairing symme- 
tries. In Eq.Q, we remove a pair from sites + r) 
and add the pair at sites (i + l,i + I + r'), whereby we 
rotate the pair by 9 degree from i ^ r direction to i — > r' 
clockwise [see Fig. ^a.)] . Note that the same results must 
be expected even if we define the spin-triplet operator as 
Ai = Ci|Ci+r| or CiiCi+ri- We also restrict ourselves to 
the case at quarter filling, i.e., n = 0.5. We calculate 
the pair-correlation functions ^ of the three-leg ladder 
with the DMRG method. The OBC are applied in the 
a;-direction so that the correlation functions will be cal- 
culated using distances taken about the midpoint. In 
the following, we study a ladder with 32 x 3 sites with 
keeping m = 4500 density-matrix eigenstates to build the 
DMRG basis. The obtained ground-state energy is ex- 
pected to be accurate to parts in 10~^t. In Fig.[3Ub), the 
pair-correlation functions D{9,1) are plotted as a func- 
tion of distance / for all possible 9 values with a fixed 
anisotropy t' — 0.5. We can see that D{9,1) decay as 
~ 1/P (1 < 7 < 2) for all 9 values. The exponent 7 
smaller than 2 would imply an attractive interaction be- 
tween electrons 0. It is consistent with our previous 
results [4^. We also find that the values of D(9,l) are 
mostly positive for 9 = 0°, 120°, 240° and, whereas, neg- 
ative for 9 = 60°, 180°, 300°. In other words, the pair 
wave function changes its sign by 7r/3 rotation. 

In order to organize the angle dependence of D{9,1) 
more explicitly, we sum it up for pair separations / = 
7-25: 



(3) 



1=7 



The open-end effects are expected to be rather small for 
central 25 sites of the system (Z ^ 25). In FiglH the 
results of D{9) with t' = 0.25 and 0.5 are shown. For 
both parameters we can see that the sign of D{9) changes 
by 7r/3 rotation, which clearly indicates the /-wave spa- 
tial symmetry of the pair correlation function. Since the 
spin-triplet pairs are formed mostly on the zigzag bond 
for t' = 0.25 [i^l, D{9) involving the on-leg pairs, namely, 
at 9 = 27r/3,57r/3 in the upper left panel, 9 = 7r/3,47r/3 
in the upper right panel, and all of 9 in the lower panels, 
are relatively small. When t' is increased from 0.25 to 0.5, 
the pairing correlation on the zigzag bond is somewhat 
reduced and, whereas, that on the leg bond is enhanced. 
The situation is rather complicated: t' enhances the an- 
tiferromagnetic correlation on the zigzag bond by the di- 
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60 120 180 240 300 60 120 180 240 300 360 
6 [degree] 

FIG. 4: (Color online) Rot at ion- angle dependence of the sum 
of the spin-triplet pair-correlation functions D{9) for J = J' — 
0. Inset shows the position of the pair annihilated by the pair 
operator Ai. 



rect exchange interaction (scaled as oc and also en- 
hances the ferromagnetic correlation on each triangle by 
the three-site cyclic-hopping exchange interaction (scaled 
as oc tt'^). Consequently, the /-wave symmetry of the su- 
perconducting state becomes more isotropic. However, 
we note that the spin-triplet pair correlation becomes 
smaller as the system approaches to an isotropic triangle 
lattice. 

In summary, the ground state of the 2D triangular- 
lattice Hubbard model is studied with the DMRG 
method. We find that the three-site cyclic-hopping mech- 
anism really induces ferromagnetic correlation around 
quarter filling, which is mandatory for the spin-triplet 
pairing, provided that the anisotropy is rather strong. 
Slow decay of the pair-correlation functions with dis- 
tance is indicative of the dominant spin-triplet supercon- 
ductivity. The pair-correlation function changes its sign 
by 7r/3 rotation, indicating the /-wave rotational sym- 
metry. Thus, we suggest a new mechanism of the /-wave 
spin-triplet superconductivity derived from the three-site 
cyclic- hopping ferromagnetic interactions. We also point 
out that the mechanism may be possibly relevant to the 
spin-triplet superconducting systems, such as Bechgaard 
salts (TMTSF)2X, cobalt oxide Nao.35CoO2-l.3H2O, and 
layered perovskite Sr2Ru04. 
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